Draft version February 1, 2008 

Preprint typeset using lAT^^X style cmulateapj v. 04/03/99 



ROSAT EVIDENCE FOR INTRINSIC OXYGEN ABSORPTION IN COOLING FLOW 

GALAXIES AND GROUPS 

David A. Buote^ 

UCO/Lick Observatory, University of California at Santa Cruz, Santa Cruz, CA 95064; buote@ucolick.org 
Accepted for Publication in The Astrophysical Journal 

ABSTRACT 

The existence of large quantities of gas that have cooled and dropped out of the hot phase in massive 
elliptical galaxies, groups, and clusters is the key prediction of the inhomogeneous cooling flow scenario. 
Using spatially resolved, deprojected ROSAT PSPC spectra of 10 of the brightest cooling flow galaxies 
and groups with low Galactic column densities we have detected intrinsic absorption over energies ~ 
0.4 — 0.8 keV at the 2a /3a level in half of the sample. Since no intrinsic absorption is indicated for 
energies below ^ 0.4 keV, the most reasonable model for the absorber is coUisionally ionized gas at 
temperatures T = 10^^® K with most of the absorption arising from ionized states of oxygen but with 
a significant contribution from carbon and nitrogen. The soft X-ray emission of this warm gas can 
also explain the sub-Galactic column densities of cold gas inferred within the central regions of most 
of the systems. (This could not be explained by an absorber composed only of dust.) Attributing the 
absorption to ionized gas reconciles the large columns of cold H and He inferred from Einstein and ASCA 
with the lack of such columns inferred from ROSAT . 

Within the central ^ 10 — 20 kpc, where the constraints are most secure, the mass of the ionized 
absorber is consistent with most (perhaps all) of the matter deposited by a cooling flow over the lifetime 
of the flow. Since the warm absorber produces no significant H or He absorption the large absorber 
masses are consistent with the negligible atomic and molecular H inferred from Hi and CO observations 
of cooling flows. It is also found that if T ^ 2 x 10^ K then the optical and FUV emission implied by the 
warm gas does not violate published constraints. An important theoretical challenge is to understand 
how the warm temperature is maintained and how the gas is supported gravitationally, and we discuss 
possible solutions to these problems that would require fundamental modification of the standard cooling 
flow scenario. Finally, we discuss how the prediction of warm ionized gas as the product of mass drop-out 
in these and other cooling flows can be verified with new Chandra and XMM observations. 

Subject headings: cooling flows - intergalactic medium ~ X-rays: galaxies 

L INTRODUCTION 

The evolution of the hot gas in the centers of massive 
elliptical galaxies, groups, and galaxy clusters has been 
most frequently interpreted in terms of the cooling flow 
paradigm (e.g., Fabian 1994). However, the characteris- 
tic radially increasing temperature profiles and centrally 
peaked X-ray surface brightness profiles usually attributed 
to cooling flows can be successfully described by only as- 
suming a two-tier structure for the gravitational potential 
(e.g., Ikebe et al 1996; Xu et al 1998); i.e., the cooler gas 
sits in the shallower potential associated with the central 
galaxy whereas the hotter gas sits in the deeper potential 
of the surrounding group or cluster. The two-tier structure 
is typically incorporated into cooling flow models (e.g., 
Thomas, Fabian, & Nulsen 1987; Brighenti & Mathews 
1998), but since cooling flows are not necessarily required 
to explain the temperatures and surface brightness profiles 
of the hot gas why should they still receive attention? 

Unlike the empirical two-tier potential model, cooling 
flows attempt to offer a nearly complete theoretical de- 
scription of the time evolution of the gas properties. A 
cooling flow describes the dynamical evolution and the X- 
ray emission of the hot gas by considering the following 
simple picture for the energy balance of a parcel of gas. 
Since the parcel of gas emits X-rays it loses energy. As the 
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parcel radiates it sinks deeper (i.e., flows inward) into the 
potential well of the system where it encounters regions 
of higher density and therefore higher pressure. Conse- 
quently, the gas parcel contracts but is then heated as 
the result of PdV work. This balance between cooling and 
heating is expected to apply over much of the cooling flow. 
Within the central regions of highest density this balance 
is broken as cooling overwhelms the heating leading to 
the key prediction of the inhomogeneous cooling flow sce- 
nario: in massive elliptical galaxies, groups, and clusters 
large quantities of gas should have cooled and dropped 
out of the flow and be distributed at least over the central 
regions of the flow. 

This prediction has inspired many searches in Hi and 
CO for cold gas at the centers of cooling flows, and all such 
attempts have either detected small gas masses or placed 
upper limits which are in embarrassing disagreement with 
the large masses expected to have been deposited in a cool- 
ing flow (e.g., Bregman, Hogg, & Roberts 1992; O'Dea et 
al 1994). If instead the mass drop-out is in the form of 
dust then current constraints on the infrared emission in 
cluster cores are not inconsistent with cooling flow models 
(e.g., Voit & Donahue 1995; Allen et al 2000b). But cool- 
ing flows are not required to explain the infrared data in 
clusters ( Lester et al, 1995 ) and individual elliptical galax- 
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ies ( [Tsai fc Mathews, 199q ). 

The case for mass drop-out received a substantial boost 
with the discovery of intrinsic soft X-ray absorption in the 
Einstein spectral data of cooling flow clusters (White et 
al 1991; Johnstone ct al 1992). The Einstein results have 
been verified with multitemperature models of the ASCA 
spectral data of clusters (Fabian et al 1994; Allen et al 
2000b), elliptical galaxies, and groups (Buote & Fabian 
1998; Buote 1999, 2000a; Allen et al 2000a). If the soft 
X-ray absorption is interpreted as cold gas then the large 
intrinsic column densities of cold H suggested by the Ein- 
stein and ASCA observations still suffer from the tremen- 
dous disagreement with the H i and CO observations noted 
above. The Einstein and ASCA results appear even more 
suspect when considering that in systems with low Galac- 
tic columns (where any intrinsic absorption should be eas- 
ier to detect) no significant excess absorption from cold 
gas is ever found with the ROSAT PSPC which should 
be more sensitive to the absorption because of its softer 
bandpass, 0.1-2.4 keV (e.g., David et al 1994; Jones et al 
1997; Briel & Henry 1996). 

We have re-examined the ROSAT PSPC data of cooling 
flows to search for evidence of intrinsic soft X-ray absorp- 
tion and in particular have allowed for the possibility that 
the absorber is not cold. Previously in Buote (2000c; here- 
after PAPER2) we have presented temperature and metal- 
licity profiles of the hot gas in 10 of the brightest cooling 
flow galaxies and groups inferred from deprojection analy- 
sis of the PSPC data. We refer the reader to that paper for 
details on the data reduction and deprojection procedure. 

In this paper we present the absorption profiles of the 
10 galaxies and groups, each of which have low Galactic 
column densities (see Table 1 of PAPER2). Partial results 
for two systems analyzed in the present paper (NGC 1399 
and 5044) also appear with results for the cluster A1795 in 
Buote (2000b; hereafter PAPERl). In §|we describe the 
models used to parameterize the soft X-ray absorption. 
We present the radial absorption profiles f or a standard 
absorber model with solar abundances in S3.1. The ef- 



dial annulus. Multiphase models are examined in §3.6 



fects of partial covering and the s ensit ivity of th e results 
to the bandpass are discussed in |3.2| and §3.3. The re- 



sults of model ing the absorption with an oxygen edge are 
presented in §3.4. In i^.t we co nsider multiphase mod- 
els such as cooling flows (§3.6.1). Evidence for em ission 
from a warm gaseous component is described in § |3.6.3 



We demonstrate the consistency of the ROSAT and ASC/ . 
absorption measurements in S^. In ^ we discuss in detail 
the implications of our absorption measurements for the 
physical state of the absorber, the cooling flow scenario, 
observations at other wavelengths, and theoretical mod- 
els. Finally, in ^ we present our conclusions and discuss 
prospects for verifying our prediction of warm ionized gas 
in cooling flows with future X-ray observations. 

2. MODELS 

2.1. Hot Plasma 

As discussed in PAPER2 we use the MEKAL plasma 
code to represent the emission from a single temperature 
component of hot gas. Because of the limited energy res- 
olution of ROSAT we initially focus on a "single phase" 
representation of the hot gas such that a single temper- 
ature component exists within each three-dimensional ra- 



2.2. Absorber 

It is standard practice to represent the soft X-ray ab- 
sorption arising from the Milky Way by material with so- 
lar abundances distributed as a foreground screen at zero 
redshift. In this standard absorption model the X-ray 
flux is diminished according to A{E) = exp(— A'hc(£')), 
where iVn is the hydrogen column density and a{E) is the 
energy-dependent photo-electric absorption cross section 
for an absorber with solar abundances. We allow A'h to 
be a free parameter in our fits to indicate any excess ab- 
sorption intrinsic to a galaxy or group and also to allow for 
any errors in the assumed Galactic value for Nn and for 
any calibration uncertainties. Note that in this standard 
model A^H is measured as a function of two-dimensional 
radius, R, on the sky. 

Intrinsic absorption is expressed more generally as, 

A{E) = /exp(-AH(7[S(l + z)]) + (1 - /), 

where z is the source redshift and / is the covering factor. 
Since the redshifts are small for the objects in our sample 
any absorption in excess of the Galactic value indicated 
by the standard model is essentially that of an intrinsic 
absorber with / = 1 placed i n fr ont of the source. We 
discuss the effects of /< 1 in §3.2. 



As discussed in PAPERl we consider oxygen absorption 
intrinsic to the galaxy or group which we represent by the 
simple parameterization of an edge, exp [—t{E/Eo)^^] for 
E > Eq, where Eq is the energy of the edge in the rest 
frame of the galaxy or group, and r is the optical depth. 
To facilitate a consistent comparison to the standard ab- 
sorber model we place this edge in front of the source, and 
thus T is also measured as a function of two-dimensional 
radius, R, on the sky. Models with / < 1 behave similarly 
to the solar-abundance absorber (§3.2). 

The photo-electric absorption cross sections used in 
this p aper are giv en by Balucihska-Church fc McCammon 
|(1992| ). Although [Arabadjis fc Bregman (1999| ) point out 
that the He cross section at 0.15 keV is in error by 13%, 
since we analyze E > 0.2 keV we find that our fits do 
not change when using the Morrison fc McCammon (1983) 
cross sections which have the correct He value. I^'urther de- 
tails on the spectral models used in the analysis are given 
in PAPER2. 

3. RADIAL ABSORPTION PROFILES 

Following PAPER2 we plot in Figures ||-^ the radial pro- 
files of the column density and oxygen edge optical depth 
obtained from the deprojection analysis according to the 
number of annuli for which useful constraints on the pa- 
rameters were obtained. This categorizes the systems es- 
sentially according to the S/N of the data. We refer the 
reader to PAPER2 for the temperature and metallicity 
profiles corresponding to these models. 

In several cases the column densities and optical depths 
could not be constrained in the outermost annuli. Ow- 
ing to the nature of the deprojection method large errors 
in the outer annuli can significantly bias the results for 
nearby inner annuli. Hence, in some systems we fixed the 
column densities to their nominal Galactic values or the 
edge optical depths to zero in the relevant outer annuli. 
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Fig. 1. — Results of the deprojection analysis for the systems where interesting constraints wore obtained in 7 annuli. (Left panels) The 
column density profiles for a standard foreground absorber with solar abundances obtained from fits over 0.2-2.2 keV (i.e., -Emin = 0-2 keV) 
are denoted by open squares for best fit and solid diamonds for Icr error ba rs: models which also have an intrinsic oxygen edge are represented 
by crosses and dashed diamonds. The Galactic hydrogen column density (Dickey &; Lockman, 1990) is shown as a dotted line. Radial units 
are arcminutes on the bottom axis and kpc on the top. (Middle panels) ine column density profiles for a standard absorber obtained from 
fits over 0.5-2.2 keV (i.e., i?min = 0.5 keV). (Right panels) Optical depths for an oxygen edge at 0.532 keV (rest frame) obtained from fits 
with i?niin = 0.2 keV. The open squares and solid diamonds are respectively the best fit and Icr errors for models where the standard absorber 
has A'^H fixed to the Galactic value; crosses and dashed diamonds refer to models with variable A'^h and thus correspond to similarly marked 
models in the left panels. 



3.1. Foreground Absorber with Solar Abundances 

We begin by examining the spectral fits using the stan- 
dard absorption model of a foreground screen (z = 0) with 
solar abundances. The left panels of Figures lllfl show 



N}i{R) obtained from spectral fits over the energy range 
0.2-2.2 keV. Our column density profiles are consistent 
with those presented in previous ROSAT studies (Forman 
et al 1993; David et al 1994; Trinchicri et al 1994; Kim 
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Fig. 2. — As Figure ^ but for systems with 5 or 6 annuli. 
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& Fabbiano 1995; Jones et al 1997; Trinchieri et al 1997; 
Buote 1999) after accounting for the different plasma codes 
and solar abundances used. 

The column densities are always within a factor of ~ 2 
of the Galactic value (N^^^), though in most cases A^h 
decreases as R decreases such that A^h < A^h small 
R. The quality of the fits for four of these systems is also 
formally poor (P < 0.01) in the central 1' bin (Table 
and for several objects in the sample the metallicities are 



very large and very inconsistent with all ASCA studies 
(see PAPER2).2 Since TVr ~ iV^'"' in the outer radii the 
observation that Nh < N^'^^ at small R indicates that, 
whatever the origin of the deficit, it must be int rinsic to 
the source. (We provide an explanation in i |3.6.3 .) 

The approximately Galactic columns arc wholly incon- 
sistent with the large excess columns inferred from mul- 
titemperature models of the spatially integrated ASCA 



^Here P represents tlie null hypothesis probability under the assumption of gaussian random errors. We discuss the suitability of this 
approximation for interpreting goodness of fit in section 3.4 of PAPER2. 
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Fig. 3. — As Figure ^ but for systems with 4 annuli. 
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spectral data of these systems (Buote & Fabian 1998; 
Buote 1999, 2000a; Allen et al 2000a). We now address 
the origin of this inconsistency. 

3.2. Effects of Partial Covering 

It has been suggested that the reason why analyses 
with the PSPC do not infer large excess column densi- 
ties for cluster cooling flows is that the standard fore- 
ground model used above systematically underestimates 
the true column intrinsic to the system (Allen & Fabian 
1996; Sarazin, Wise, & Markevitch 1998). However, in 



PAPERl we have tested this hypothesis for NGC 1399 
and 5044 (and the cluster A1795) using our deprojection 
code. We find that the hot gas within the central r = 1' 
(3D) cannot be absorbed very differently from the gas pro- 
jected from larger radii because their spectral shapes for 
energies below ~ 0.5 keV are very similar. If we do assume 
an absorber with covering factor / — 0.5 we obtain an ex- 
cess column AA'^H = at best fit and AA'^h < ^h''' at 
> 90% confidence for NGC 1399 and 5044. (Note that the 
/ = 0.5 model implies a flat absorbing screen that bisects 
the source so that the 2D and 3D radii are equal, and thus 
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Fig. 4. — As Figure pJ but for the one system with 3 annuli. 



Table 1 

Quality of Fits and Edge Optical Depth for r = 1' 









No Edf 


;e 






Edge 




F Test 


Name 


Nil 




dof 


P 




dof 


P 


T 


P{F) 


N507 


Fix 


36.2 


47 


0.87 


32.2 


46 


0.94 


1.9(0.59) 


2.7E-2 




Free 


34.0 


46 


0.91 


30.4 


45 


0.95 


1.5(0.49) 


2.6E-2 


N533 


Fix 


29.9 


31 


0.52 


29.7 


30 


0.48 


0.55(0.00) 


0.66 




Free 


28.9 


30 


0.52 


28.9 


29 


0.47 


0.00(0.00) 


1.00 


N1399 


Fix 


106.9 


97 


0.23 


83.7 


96 


0.81 


3.8(0.89) 


1.3E-6 




Free 


104.9 


96 


0.25 


83.7 


95 


0.79 


3.3(0.97) 


3.9E-6 


N2563 


Fix 


24.0 


25 


0.52 


22.9 


24 


0.53 


1.5(0.00) 


0.29 




Free 


23.6 


24 


0.48 


22.9 


23 


0.47 


1.3(0.00) 


0.41 


N4472 


Fix 


209.0 


104 


4.7E-9 


158.1 


103 


4.0E-4 


1.5(0.74) 


8.8E-8 




Free 


158.9 


103 


3.4E-4 


142.2 


102 


5.3E-3 


1.3(0.25) 


7.9E-4 


N4636 


Fix 


109.9 


52 


4.9E-6 


105.8 


51 


l.OE-5 


1.3(0.00) 


0.17 




Free 


104.7 


51 


1.4E-5 


103.1 


50 


1.5E-5 


0.58(0.00) 


0.38 


N4649 


Fix 


67.2 


61 


0.27 


58.1 


60 


0.55 


1.7(0.44) 


3.2E-3 




Free 


63.4 


60 


0.36 


57.8 


59 


0.52 


1.7(0.00) 


2.0E-2 


N5044 


Fix 


191.2 


118 


2.3E-5 


156.5 


117 


8.7E-3 


1.8(0.59) 


9.2E-10 




Free 


138.6 


117 


8.5E-2 


126.4 


116 


0.24 


1.3(0.45) 


l.lE-3 


N5846 


Fix 


84.0 


47 


7.4E-4 


78.4 


46 


2.1E-3 


0.83(0.00) 


7.6E-2 




Free 


74.6 


46 


4.8E-3 


72.6 


45 


5.7E-3 


0.57(0.00) 


0.27 


H62 


Fix 


37.8 


46 


0.80 


37.7 


54 


0.77 


0.16(0.00) 


0.73 




Free 


38.6 


45 


0.74 


37.7 


44 


0.74 


0.46(0.00) 


0.31 



Note. — Results of deprojection analysis in the central radial bin {R = 1'). Models where the column density of the standard foreground 
absorber with solar abundances is fixed to the Galactic value are indicated by "Fix" under the Nn column; variable column models are indi- 
cated by "Free". For the models with an intrinsic oxygen edge at 0.532 keV (rest frame) we list the best-fitting optical depth, r, with its 95% 
confidence lower limit in parentheses. The null hypothesis probability is listed in the column P under the assumption of gaussian random 
errors (see section 3.4 of PAPER2). The F-Test probability, P(F) (e.g., Bevington 1969), quantifies the improvement in when adding the 
edge; i.e., P{F) <C 1 indicates significant improvement. 



the values of AA^h quoted do refer to quantities within the 
3D radius r — 1'.) 

Entirely analogous results are obtained for the other sys- 
tems in our sample. We mention that partial covering 
models never improve the fits over the / = 1 case. The 
only effect is that somewhat larger columns are generally 
allowed; e.g., for / — 0.5 the implied upper limits for the 
excess columns are typically a factor of 20%-40% larger 
than for / = 1. 

Hence, as in PAPERl we conclude that models with 
/ < 1 cannot account for (1) the large excess columns in- 
ferred from ASCA , (2) the sub-Galactic columns and poor 



fits obtained for several systems in the central 1', or as we 
now discuss (3) the sensitivity of A^h to the lower energy 
boundary of the bandpass. 

3.3. Sensitivity of N}i{R) to Bandpass 

Thus far we have shown (and confirmed previous results) 
that the ROSAT PSPC does not indicate the presence of 
excess absorbing material arising from cold gas intrinsic 
to the galaxy or group. In Figure |^ (Left) we display the 
ROSAT PSPC spectrum of NGC 1399 within the central 
arcminute along with the best-fitting plasma model modi- 
fied by a standard absorber with Nh = N^^^. The simple 
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Fig. 5. — (Left) ROSAT PSPC spectrum within i? = 1' of NGC 1399. Also shown is the best-fitting single-phase plasma model modified 
by the standard absorber with Nn = N^'^^. (Right) The best fitting model with variable A^h when only energies above 0.5 keV are included 
in the fit. In this case the fitted column density is A'^h ~ 15Af^^'; i.e., when Emin = 0-5 keV the standard absorber model predicts that there 
should be negligible emission at lower energies in conflict with the observation. 



model provides a good visual fit to these data as well as a 
formally acceptable fit (x^ = 107 for 97 dof and P = 0.23). 

However, the evidence for intrinsic soft X-ray absorp- 
tion from cold gas in galaxies, groups, and clusters from 
the Einstein SSS and the ASCA SIS is obtained with data 
restricted to energies above ~ 0.5 keV because lower en- 
ergies reside outside the bandpasses of those instruments. 
Let us examine what happens to the spectral fits of the 
PSPC data of NGC 1399 within the central arcminute if 
instead we raise the lower energy limit of the bandpass, 
£^min, to a value near 0.5 keV comparable to ASCA and 
Einstein , and we allow TVjj to be a free parameter. 

We find that when i^min ~ 0.2—0.3 keV we obtain fits es- 
sentially as indicated in Figure || (Left) with iVn « N^'^^. 
We see a noticeable change near Smin — 0.4 keV when 
the fitted column density increases to iVn ~ 2iV^^'. A 
dramatic increase occurs near -Bmin — 0.5 keV which 
we show in Figure ^ (Right). The best-fitting model 
gives Nu « 157V^''i when data below 0.5 keV are ex- 
cluded from the fits: the standard absorber model when 
£^min = 0.5 keV predicts that there should be negligible 
emission at lower energies in clear conflict with the data 
for E ^ 0.2 — 0.3 keV. The large column density obtained 
when i^niin = 0.5 keV is similar to the large value inferred 
from ASCA with two-temperature or cooling flow models 
(see §5 in Buote 1999). 

For larger values of Emm we flnd that A^h does not 
change significantly within the uncertainties. The statisti- 
cal uncertainties on the fitted spectral parameters increase 
with increasing Emin since the degrees of freedom decrease 
as the lower energy data are ignored. 

In the middle panels of Figures |l]^ we plot iVH(i?) for 
^-min — 0.5 keV. The character of the Ah profiles for 
^-min = 0.5 keV is entirely different from the previous 
Emin = 0.2 keV case for half of the sample: NGC 507, 
1399, 4472, 4649, and 5044. In these systems AH(i?) for 
^-min — 0.5 is consistent with the Galactic values in the 
outermost annuli and increases as R decreases until (ex- 
cept for NGC 4472) it reaches a value consistent with a 
maximum for i? ^ 1'. For the other five galaxies the con- 



straints are too poor to discern a trend, though within the 
large errors their profiles are consistent with increasing as 

The excess column densities inferred when Emin = 0.5 
keV are most significant for NGC 1399. The lower lim- 
its within i? = 1' are 8 X 10^° cm^^ and 3 x 10^° cm'^ 
at 95% and 99% confidence respectively which are factors 
of « 6 and 3 larger than the Galactic value. Also within 
i? = 1' for NGC 507 we obtain 95%/99% lower limits of 
8/6 X 10^" cm^^ compared to the adopted Galactic value 
of 5.2 X 10^*^ cm^^. The most significant measurement for 
NGC 5044 is within the i? = 1' - 2' annulus where we find 
that A^H > Af^""' at the 95% confidence level. 

In order to obtain measurements of the intrinsic absorp- 
tion at a higher significance level we must include the data 
below 0.5 keV. This requires a more appropriate model of 
the intrinsic absorption that does not confiict with the 
emission at lower energies which we now consider. 

3.4. Intrinsic Oxygen Edge 

Since we find that A'H(-Emin) ~ constant for -Emin ^ 0.5 
keV, the portion of the spectrum responsible for the excess 
absorption must be near 0.5 keV. Considering the PSPC 
resolution [AE/E = 0.43(£;/0.93keV)-"-^] and effective 
area this translates to energies ~ 0.4 — 0.7 keV. The dom- 
inant spectral features in both absorption and emission 
over this energy range are due to oxygen, thoug h io nized 
carbon and nitrogen can contribute as well (see § |5.2D . 

Since the PSPC data cannot distinguish between a sin- 
gle edge and multiple edges, we parameterize the intrin- 
sic absorption with a single absorption edge at 0.532 keV 
(rest frame) corresponding to cold atomic oxygen (Oi). In 
the right columns of Figures we plot the optical depth 
profiles, t{R), for the Oi edge obtained from fits with 
Emin = 0.2 keV. 

For every system we find that the shape of t{R) is very 
similar to that of A^h(-R) for Emin = 0.5 for the standard 
absorber. Moreover, for those systems where A'h ^ A^^' 
for Emin — 0.5 we find that t{R) ^ 0.1 in the outer- 
most annuli and increases to t(R) ~ 1 in the central bin. 
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Therefore, the single oxygen edge reproduces all of the ex- 
cess absorption indicated by Nii{R) for ii^niin = 0.5 for the 
standard absorber model. 

The fits within the central radial bin are clearly im- 
proved for several systems when the single oxygen edge is 
added, even though in some cases the quality of the fit is 
already judged to be formally acceptable (null hypothesis 
P ^ 0.1) without the edge. The improvement in the 
fit is quantified by the F Test (e.g., Bevington 1969), and 
in Table ^ for each galaxy and group we give the F-Test 
probability, P{F), which compares the fits with and with- 
out the edge; i.e., P{F) ^ 1 indicates the edge improves 
the fit significantly. 

Concentrating on models with the standard absorber 
with Nil — N§^^ we see in Table |^ the largest improve- 
ments exist for NGC 5044 (P(F) = 9.2 xlO^^"), NGC 4472 
{P{F) = 8.8 X 10"*^), and NGC 1399 (P(F) = 1.3 x 10"^). 
A very significant improvement is also found for NGC 4649 
{P{F) = 3.2 X 10"'^) while marginal improvements are 
found for NGC 507 (P(F) = 2.7 x 10"^) and NGC 5846 
[P{F) = 7.6 X 10-2). Ahhough only NGC 4472, 5044, 
and 5846 are indicated to have formally unacceptable fits 
in terms of the '^uU hypothesis probability for models 
with or without the edge, the fact that adding the edge 
lowers by much more than 1 (i.e., P{F) <C 1) in several 
cases indicates that the fits are indeed better with the edge 
in those systems. 

For every system whe re we found Nn < N^^^ for the 
standard absorber (§3T), we find that Nn systematically 
increases when the oxygen edge is added. Although the 
addition of the edge results in A^h ~ N^^^ for many of 
these systems we find that in some cases (most notably 
NGC 5044) A^H is still significantly less than A^^^\ and 
the fits in the central bin, though improved, are still for- 
mally unacceptable. Thus, adding the edge does signifi- 
cantly improve the models, but it apparently is not the 
only improvement required in some cases. We discus s an- 



other mechanism to improve the fits below in §3.6.3 



We obtain the best constraints on the edge optical depth 
for NGC 507, 1399, 4472, 4649, and 5044. Only for these 
systems is r > significant at the 95% confidence level in 
the central bin whether or not A^h is fixed to the Galac- 
tic value for the standard absorber model. In some cases 
r > is significant at the 99% level - we discuss individual 
cases below in 
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We have also mvestigated whether these oxygen absorp- 
tion profiles can be reproduced by a profile of decreasing 
oxygen abundance; i.e., the strong Ka lines of O VII and 
O VIII lie between 0.5-0.65 keV. If instead we allow the 
oxygen abundance in the hot gas to be a free parameter 
in the fits we find that typically the best-fitting oxygen 
abundance is zero, and the quality of most of the fits are 
improved similarly to that found when the oxygen edge is 
added. This degeneracy is not surprising owing to the lim- 
ited energy resolution of the PSPC. The notable exception 
is NGC 1399 where the fits for r — V are only improved 
to = 95.9 for zero oxygen abundance as opposed to 
83.7 for the edge (both have variable A'h). However, zero 
oxygen abundance in the centers is highly unlikely because 
of the expected enrichment from the stars in the central 
galaxy. And if instead we consider plausible 0/Fe ratios 
to be at least 1/2 solar in NGC 507, 1399, 4472, 4649, 



and 5044, then the fits are not as much improved as when 
adding the edge. 

As mentioned above (and in PAPERl) we find that for 
most radii the constraints on the edge energy are not very 
precise which is why we fixed the edge energy in our anal- 
ysis. The best constraints are available for NGC 1399 and 
5044 within the central bin for which we obtain 0.5llg'Q5 

keV and 0.5ltJ5;|5^ keV (90% confidence) for the edge en- 
ergies for the standard absorber models with variable A^h. 
(Models with fixed foreground Galactic columns give sim- 
ilar resuhs; e.g., 0.53l[5;5]^ keV for NGC 1399.) 

These constraints are consistent with the lower ioniza- 
tion states of oxygen but not edges from the highest states 
Ovi-viii. Due to the limited resolution we can add addi- 
tional edges to "share" the r obtained for the Oi edges, 
although when using a two-edge model a significant t can- 
not be obtained for edge energies above ~ 0.65 keV corre- 
sponding to ~Ovi. 

3.5. Comments on Individual Systems 

We elaborate further on the results for individual sys- 
tems. When comparing A^h profiles to those obtained from 
previous ROSAT PSPC studies we implicitly account for 
any differences in the plasma codes and solar abundances 
used. 

3.5.1. Systems with 7 Annuli 

In Figure |^ we display the results for the three systems 
where the spectral parameters are well determined in seven 
annuli. These observations thus generally correspond to 
the highest S/N data in our sample. The evidence for 
intrinsic oxygen absorption is strongest for these systems. 

NGC 507: The oxygen edge optical depth, r(i?), falls 
gradually from the center and remains significantly non- 
zero out to i? ~ 4'; e.g., in the R = 3'-4.25' annulus the 
95% confidence lower limits on r are 0.39 and 0.24 respec- 
tively in models with fixed (Galactic) and variable A^h for 
the standard absorber. The values of A'h for the stan- 
dard absorber are consistent with those obtained by Kim 
fc Fabbiano (1995] ) with the PSPC data. 
NGC 1399: This system exhibits the most centrally 
peaked r profile in our sample. Within the central ar- 
cminute r > 0.26 and 0.33 at 99% confidence respectively 
for the fixed and free A^h models. At larger radii the 
non-zero optical depths are also quite significant; e.g., for 
R = 2.5'-4' T > 0.64 and 0.50 at 99% confidence for the 
fixed and free A^h models. We obtain values of A'h for 
the stand ard absorber cons istent with the previous PSPC 
study by Jones et al (1997). 

NGC 5044: The second radial bin {R = l'-2') actually 
has a smaller uncertainty on the oxygen edge optical depth 
than the central bin; i.e., the 95% lower limits on t are 
0.87 and 0.79 respectively for the fixed and free A^h mod- 
els. In fact, the corresponding 99% lower limits are 0.75 
and 0.42 for R = l'-2' which are larger than the values for 
the R — 1' bin. The optical depth remains significantly 
non-zero out to the R = 3'-4.5' bin in which we obtain 
95% lower limits on r of 0.36 and 0.10 for the fixed and 
free A^h models. Our values of A^h are consistent wit h 
those obtained from the PSPC data by David et al (1994). 



3.5.2. Systems with 5-6 Annuli 
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In Figure Q we display the results for the three systems 
where the spectral parameters are well determined in 5-6 
annuli. Only for NGC 4472 is the intrinsic oxygen absorp- 
tion clearly significant. 

NGC 2563: The uncertainties on t{R) are large and con- 
sistent with zero at the center. However, the shapes of the 
error regions, especially the large uncertainty at the center, 
are consistent with the same type of increasing r profile 
with decreasing R found for the systems with seven annuli. 
NGC 4472: The optical depth is consistent with t ^ 2 for 
R ^ 5' and then decreases rapidly at larger radii. Unlike 
the other systems with evidence for intrinsic absorption 
r is most significant away from the central bin; i.e., for 
R = 3.25'-5' the 95% confidence lower limits on r are 1.62 
and 1.75 respectively for the fixed and free Nn models. 
(The corresponding 99% lower limits are 1.27 and 1.46.) 
Inter estingly, R 5 ' corresponds to the region where fr\ 



win & Sarazin (1996) identified holes in the X-ray emission 
from visual examination of the PSPC image, and thus the 
large oxygen edge optical depths could be related to these 
holes. The sub-Galactic columns obtained for the stan- 
dard absorb er with variable TVh are consistent with those 
obtained by Forman et al (1993). 

NGC 5846: The uncertainties on t{R) are large, and 
although a substantial amount of intrinsic absorption is 
allowed by the data, no excess absorption is required. 

3.5.3. Systems with 4 Annuli 

In Figure |^ we display the results for the three systems 
where the spectral parameters are well determined in 4 an- 
nuli. The uncertainties on t{R) are large for each of these 
systems and thus no intrinsic oxygen absorption is clearly 
required by the data. 

3.5.4. Systems with 3 Annuli 

In Figure |4| we display the results for the one system 
where the spectral parameters are well determined in only 
3 annuli, NGC 4649. In contrast to the systems with 4 
annuli we find evidence for significant intrinsic oxygen ab- 
sorption in the central bin for NGC 4649. 

NGC 4649: In the central radial bin the 90% lower hmits 
on T are 0.55 and 0.20 respectively for for the fixed and 
free TVh models, although only the 95% lower limit for 
the model with fixed A'h is significantly larger than zero 
(Table p. Our values of iVn for the standard absorber 
are consistent with thos e obtained from the PSPC data 
by [Trinchieri et al (199"^ ). 



3.6. Multiphase Models 

3.6.1. Simple Cooling Flow 

In the inhoniogeneous cooling flow scenario the hot gas is 
expected to emit over a continuous range of temperatures 
in regions where the cooling time is less than the age of 
the system. We consider a simple model of a cooling flow 
where the hot gas cools at constant pressure from sonic up- 
per temperature, Tmax (e.g., Johnstone et al 1992). The 
differential emission measure of the cooling gas is propor- 
tional to M/A(T), where M is the mass deposition rate of 
gas cooling out of the flow, and A(T) is the cooling func- 
tion of the gas (in our case, the MEKAL plasma code). 



Since the gas is assumed to be cooling from some upper 
temperature Tmax, the cooling flow model requires that 
there be a reservoir of hot gas emitting at temperature 
Tinax but is not cooling out of the flow. Consequently, our 
cooling flow model actually consists of two components, 
CF-l-lT, where "CF" is the emission from the cooling gas 
and "IT" is emission from the hot ambient gas. We set 
Tmax of the CF component equal to the temperature of 
the IT component, and both components are modified by 
the same photoelectric absorption. This simple model of 
a cooling flow is appropriate for the low energy resolution 
of the ROSAT PSPC and has the advantage of being well 
studied, relatively easy to compute, and a good fit to the 
ASCA data of many elliptical galaxies and groups (e.g., 
Buote & Fabian 1998; Buote 1999, 2000a). 

When fitting this cooling flow model to the ROSAT 
spectra of the galaxies and groups in our sample we find 
that if only absorption from the standard cold absorber 
model with solar abundances is included then we obtain 
results identical to the single-phase models; i.e., the CF 
component is clearly suppressed by the fits. Since the CF 
model includes temperature components below 0.5 keV 
it has stronger O VII and O VIII lines than the single- 
phase models which, as shown above, already predict too 
much emission from these lines. It is thus not surprising 
that only when an intrinsic oxygen absorption edge is in- 
cluded in the fits can we obtain a significant contribution 
from a CF component. 

Even when adding the oxygen edge the cooling flow 
models do not improve the fits perceptively in any case. 
And since the magnitude of t for the oxygen edge is 
degenerate with M of the cooling flow component, the 
constraints on both parameters are quite uncertain. For 
NGC 5044, which has the best constraints, adding the 
CF component improves from 156.5 to 154.1 for 116 
dof within the central arcminute (i.e., marginal improve- 
ment). The oxygen edge optical depths are consistent with 
those obtained from the single-phase analysis. Only within 
r = 2' is there an indication of significant cooling with 
M ^ 12A f0yr~^ which is ve ry consistent with the ROSAT 
results of David et al (1994 ) within the same radius. (Our 
mass deposition rates have large statistical uncertainties 
because of the inclusion of the intrinsic oxygen edge.) 

Thus, simple cooling flow models give results that are 
entirely consistent with the single-phase models. 



3.6.2. Two Hot Phases 

A two-temperature model (2T) is a more flexible multi- 
phase emission model than the constant-pressure cooling 
flow and can very accurately mimic a cooling flow spec- 
trum over typical X-ray energies (Buote, Canizares, & 
Fabian 1999). If we restrict the temperatures of the 2T 
model to lie between ~ 0.5 — 2 keV appropriate for hot gas 
near the virial temperatures of these galaxies and groups, 
then we obtain results equivalent to the cooling flow mod- 
els above; e.g., (1) the extra temperature component is 
only allowed if the intrinsic oxygen edge included, and (2) 
results are entirely consistent with the single-phase mod- 
els. Since the temperatures of the 2T models are fitted 
separately, constraints on the 2T models are even poorer 
than the cooling flows because of the extra free parameter. 
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3.6.3. Two-Phase Medium: Warm and Hot Gas 



Recall that when fitted with only a standard absorber 
with solar abundances that Nh{R) decreases as R de- 
creases such that A^H is less then N^^^ near i? = for 
most of the galaxies and groups (§^). This trend implies 
the existence of excess soft X-ray emission above that pro- 
duced by the hot gas, the signature of which is most evi- 
dent at the centers of these galaxies and groups. Although 
the sub-Galactic values of N-^ are only marginally signif- 
icant in many cases, they are highly significant for NGC 
4472 and NGC 5044. 

If the excess soft X-rays near the centers represent emis- 
sion from coronal gas, then the temperatures must be 
< 0.1 keV (i.e., distinctly cooler than the hot gas phase). 
For NGC 5044, which has the most significant sub-Galactic 
column densities in the central bins, we find that when 
adding an extra temperature component to the single- 
phase model modified only by the standard absorber with 
TVh = iV^^' (i.e., no intrinsic oxygen edge), that the fits 
are improved very similarly to the single-phase models 
when A^H is allowed to take a value significantly less than 
N^'^^. For example, in the central bin the fit improves 
from = 191 (118 dof) for the IT model to = 133 
(116 dof) for the 2T model (both models with iVn = N^""^) 
very similar to the result obtained for the IT model with 
iVn < N^""^ (i.e., = 139 for 117 dof - "Free" in Ta- 
ble |l ) . Essentially the same large improvement is found 
for the second radial bin (i? = 1' — 2'). Significant but 
smaller improvement in the fits (Ax^ ^ 10) is also seen in 
in the third and fourth bins (i.e. R — 2' — 4.5'). 

The inferred temperature of the soft component in the 
central bin is T = 0.05^^;^^ keV {6^1 x 10^ K) at 90% 
confidence consistent with values obtained for other radii; 
i.e., consistent with "warm" gas rather than "hot" gas near 
the virial temperatures of the galaxies and groups. Gas at 
these warm temperatures is not optically thin to photons 
with energies near 0.5 keV (e.g., Krolik & Kallman 1984), 
and thus this warm gas may be res pon sible for the in- 
trinsic oxygen absorption inferred in §3.4. However, since 
this warm gas is apparently partially photoionized by the 
hot gas (and perhaps by itself) a proper calculation of the 
emission from the warm gas must consider radiative trans- 
fer effects which is beyond the scope of this paper. (We 
shall continue to use the optically thin models in this pa- 
per.) We discuss the properties of this warm gas in more 
detail in §||. 

The improvement obtained in the fit when adding a 
component of warm gas to the single-phase model (iVn = 
iV§^') is somewhat larger than that obtained when adding 
the single oxygen edge (i.e., — 133 for 116 dof for the 
2T model versus x^ = 156.5 for 117 dof for the edge - Ta- 
ble H) . Unfortunately, due to the limited energy resolution 
of the PSPC it is very difficult to obtain simultaneous con- 
straints on both absorption and emission models for both 
the warm and hot gas. 

The most reliable constraints on the warm gas compo- 
nent are possible for NGC 5044 because for IT models we 
find that TVh of the standard absorber is well below iV^^' 
at small R whether or not an intrinsic oxygen edge is in- 
cluded (Figure |l]). For a 2T model with an intrinsic oxygen 
edge we find that the temperature of the warm component 

^See bttp:/ /hcasarc. gsfc.nasa.gov/docs/rosat 



in the central bin is T = O.OGIJS;^^ keV (7^5 x lO""^ K) at 
90% confidence consistent with that obtained for the 2T 
model without the oxygen edge. The oxygen edge opti- 
cal depth is not very well constrained, t = 0.8lQ g, 1.21q'7 
(90% confidence) in the inner two bins respectively, which 
is about half the best-fitting value and near the 95% lower 
limit of T obtained without including the emission from the 
warm gas (Table |l]). The weakened constraint on the oxy- 
gen edge optical depth for the 2T model of NGC 5044 re- 
fiects the relatively small improvement in the fit (x^ of 127 
for 115 dof) over the 2T model without an edge. For the 
other galaxies, most notably NGC 1399, the constraints 
on the oxygen edge are not weakened nearly as much. 

Emission from such a warm gas component is consistent 
with, though not as clearly required by, the PSPC data 
for the other systems. Only for NGC 4472 is clearly sig- 
nificant improvement found when adding both the warm 
gas component and the oxygen edge. Temperatures (T ~ 
(5 — 10) X 10^ K) and the reduction in edge optical depth 
are obtained similarly as found for NGC 5044. Future ob- 
servations with better energy resolution are required to 
definitively confirm and measure the emission and absorp- 
tion properties of the warm gas in all of these systems. 

3.7. Caveats 

(i) Calibration: The gain of the PSPC is well calibrated 
and in particular no significant calibration problems near 
0.5 keV have been reported.'^ The large values of r ~ 1 
obtained in the central regions imply that the absorption 
significantly affects a large energy range comparable to the 
energy resolution of the PSPC; e.g., the Oi edge absorbs 
25% of the flux at 0.8 keV for r = 1; i.e. possible residual 
calibration errors near 0.5 keV where the effective area is 
changing rapidly (e.g.. Figure 1 of Snowden et al 1994) 
can not explain the need for absorption at higher ener- 
gies. Moreover, the shape of r(i?) is not the same for each 
system as would be expected if calibration were responsi- 
ble for the intrinsic oxygen absorption found in half of the 
sample; e.g., r(i?) for NGC 1399 is much more centrally 
peaked than for NGC 4472 or the other systems. The 
agreement between the absorption inferred by the PSPC 
and ASCA mentioned below in §^ further argues against 
a systematic error intrinsic to the PSPC being responsible 
for the measured oxygen absorption. 
(a) Galactic Columns: All of the objects in our sample 
reside at high Galactic latitude, and thus the Galactic ab- 
sorption should be fairly uniform over the relev ant 5'-10' 
scales. Errors in the assumed Galactic columns ( Dickey ^ 
[Lockman, 199C| ) should only affect the baseline value and 
not the variation with radius. 

(Hi) Background: Errors in the background level affect 
most seriously the lowest energies (^0.4 keV) which are 
also most sensitive to the column density of the stan- 
dard absorber model. Hence, measurements of iVu at 
the largest radii (which have lowest S/N) are the most 
affected by background errors. The hydrogen column den- 
sities measured in the outer radii (see Figures |l|-||) are very 
similar to and are usually consistent with the assumed 
Galactic values within the estimated la errors which at- 
tests to the accuracy of our background estimates. The 
intrinsic oxygen optical depths in the central radial bins 
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are very insensitive to the background level. 

4. COMPARISON TO ASCA 

The intrinsic oxygen absorption indicated by the PSPC 
data in half of our sample is most significant within the 
central 1' — 2' which is similar to the width of the ASCA 
PSF. In addition, since the ASCA SIS is limited to £; > 0.5 
keV, and the efficiency near 0.5 keV is severely limited 
due to instrumental oxygen absorption, it cannot be ex- 
pected that ASCA can distinguish an oxygen edge from a 
standard absorber with solar abundances. However, it is 
instructive to examine the consistency between results ob- 
tained from spatially resolved (PSPC) and single-aperture 
{ASCA ) methods. As mentioned in PAPERl the results 
obtained when adding an oxygen edge to the ASCA data 
of NGC 1399 and 5044 are consistent with those obtained 
with the PSPC. Since similar results are found for other 
objects in our sample showing intrinsic absorption, we fo- 
cus on NGC 1399 for illustration. 

Previously we (Buote, 1999) have fitted two- 



tempcrature models to the accumulated ASCA SIS and 
CIS data within i? « 5' of NGC 1399 and obtained 
Arc ^ 49+6 ^ -^q20 (,j^^-2 (9q% confidence) for the standard 

absorber on the cooler temperature component. Using 
the meteoritic solar abundances (see §4.1.3 of PAPER2) 
slightly modifies this result to = 401^ x 10^" cni"^ 
which is about 30 times larger than the Galactic value 
{N^ — 1.3 X 10^° cm~^). Comparing this result to those 
obtained from the PSPC for E^^in = 0.5 keV (Fig|) we see 
that the ASCA column density is qualitatively similar to 
the value of Nn « 20 x 10^° obtained within R ~ 1' and is 
consistent with the total column within R « 5'. If instead 
the columns of the standard absorber are fixed to Galactic 
on both components, and an intrinsic Oi edge is added to 
the cooler component, then we obtain r = G.Olgy (90% 
confidence) for the ASCA data. Again this ASCA result 
is similar to r 4 in the central PSPC bin and is very 
consistent with the value of r = 5.7 obtained from adding 
up the best-fitting values obtained with the PSPC within 
R = 5'. Therefore, the oxygen edge provides as good or 
better description of the excess absorption inferred from 
multitemperature models of ASCA data within the cen- 
tral few arcminutes as a standard cold absorber with solar 
abundances, and yields optical depths that are consistent 
with those obtained with the PSPC data. 

Although a similar consistency is achieved for NGC 4472 
the interpretation of the absorption of the two-component 
model of the single-aperture ASCA data is not so straight- 
forward because the absorption is not obviously concen- 
trated at the center (i.e. on the cooler temperature compo- 
nent). In fact, a spatially uniform absorber within i? = 5' 
is probably a better description of the PSPC data. If the 
column densities on both temperature components are tied 
together for the two-temperature model of the ASCA data 
the quality of the fit is slightly worse (Ax^ = 8 for 361 
dof), but A^H « 11 X 10^*^ cm^^ which is a fair representa- 
tion of the average A^h profile obtained from the PSPC for 
Emhi = 0.5 keV (Fig||). Similar agreement is found for the 
O I edge when applied to both temperature components. 



5. WARM IONIZED GAS IN COOLING FLOWS 



We consider now in some detail the properties of the in- 
trinsic absorber and their implications. Initially we focus 
our attention on the physical state of the absorber. 

5.1. Why Not Dust? 

Models of dust grains indicate that dust can give rise to 
significant amounts of absorption between 0.1-1 keV (e.g., 
Laor & Draine 1993). In principle such grains could ex- 
plain the intrinsic absorption we have inferred for energies 
above 0.5 keV in the PSPC data of half of the galaxies 
and groups in our study as well as the excess absorption 
detected for energies above ~ 0.5 keV in the ASCA data 
of bright galaxies, groups, and clusters. 

However, dust should also heavily absorb X-rays with 
energies between 0.1 — 0.4 keV (e.g., Laor & Draine 1993) 
which is inconsistent with the ROSAT data for the 10 
galaxies and groups in our stud y and, e.g., the large sa mple 
of ROSAT clusters studied by jAUen fc Fabian (1997| ). To 
evade this constraint (i.e., no absorption below 0.4 keV) 
an unconventional model for the formation of dust grains 
has to be postulated; e.g., dust grains conde nse out of 
a medium in wh ich helium remains ionized (Arnaud & 
Mushotzky, 1998). Even if we consider the (rather un- 
likcly) possibility that dust does not absorb X-rays below 
~ 0.4 keV, we still have that dust cannot account for the 
excess X-ray emission at tho se e nergi es imp lied by the sub- 
Galactic column densities (§3.1 and 3.6.3). 

Other strong arguments against the existence of large 
quantities of oxygen-absorbing dust in cooling flows have 
been made in the past. At the centers of cooling flows 
where the gas density is largest and the soft X-ray ab- 
sorption is most significant dust should not be present in 
large quantities because the grains arc rapidly destroyed 
by sputtering by the hot gas (e.g., Tsai & Mathews 1995; 
Voit fc Donahue 1 995). A different argument due to Voit 
fc Donahue (1995| ) considers that the transient heating of 
the grains by X-rays from the hot gas prevents CO from 
fully condensing onto dust grains. Consequently, signif- 
icant amounts of oxygen would remain in molecular gas 
which would produce CO emission from rotational lines 
that are inconsistent with the generally negligible CO de- 
tections in cooling flows (e.g., Bregman, Hogg, fc Roberts 
1992; O'Dea et al 1994). 

5.2. Temperature of the Warm Ionized Gas 

Unlike dust all of the key features of the observed soft 
X-ray absorption and emission can be accounted for if 
the absorber is (primarily) coUisionally ionized gas. The 
lack of intrinsic absorption observed for energies between 
- 0.2 - 0.4 keV in the ROSAT data requires that H 
and He be completely ionized implying a gas temperature 
T ^ 1.0 X 10^ K (e.g., Sutherland fc Dopita 1993). In order 
to have significant absorption at 0.5 keV the temperature 
cannot be larger than « 1 x 10^ K (see Figure 2 of Kro- 
lik fc Kallman 1984). This absorbing temperature range 
{T = 10^~^ K) is entirely consistent with that inferred 
from the gas emission to explain the sub-G alactic column 
densities especially for NGC 5044 (§ 3.6.3| ). This consis- 
tency of absorption and emission properties lends strong 
support to the idea that the absorber is warm ionized gas. 

At these temperatures carbon and nitrogen are not com- 
pletely ionized (e.g., Nahar fc Pradhan 1997; Sutherland 
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& Dopita 1993), and thus these elements wiU also con- 
tribute to the soft X-ray absorption. For nitrogen the 
states Niv-vi are significant, and their edge energies span 
0.46-0.55 keV (rest frame) which are consistent with the 
edge energy range determined from the PSPC data. Since 
N/0 = 1/8.51 (assuming solar abundance ratios), and the 
threshold cross sections for absorption are similar for N 
and O, only ~ 12% of the optical depth we have measured 
in each system likely arises from ionized nitrogen. 

The ionization fraction for carbon changes rapidly near 
T = 10^ K with Cv dominating for temperatures above 
this value and up to T « 10^ K. The edge energy for Cv is 
0.39 keV, and since C/0 « 0.5 (assuming solar abundance 
ratios) and the threshold cross sections of C and O are sim- 
ilar, the optical depth of C v is about half that expected 
from a dominant ionized state of oxygen. However, the 
strong instrumental carbon absorption leaves the PSPC 
with essentially no effective area over energies 0.28 keV to 

0.4 keV (e.g., see Figure 1 of Snowden et al 1994), and 
thus it is only possible to detect intrinsic absorption from 
the Cv edge for energies above ~ 0.4 keV even consid- 
ering the smearing due to the limited energy resolution. 
This is entirely consistent with the variation of A'h with 



-Bmin described in 53.3 



Considering the energy resolution and effective area 
curve of the PSPC, the 0.532 keV edge that we have used 
to parameterize the intrinsic absorption is likely an average 
of the ~ 40% contribution from the Cv edge (0.39 keV) 
and Niv-vi edges (0.46-0.55 keV) with a ~ 60% contribu- 
tion from ionized oxygen states. Although as discussed in 
§ ^.4| the PSPC data cannot distinguish between multi-edge 
models, when using a more realistic absorber model con- 
sisting of C V and N vi and an oxygen edge, we are able to 
obtain a significant optical depth for the oxygen edge for 
energies as high as ~ 0.75 keV consistent with Ovii (0.74 
keV). However, to insure that oxygen produces at least as 
much absorption as the C and N edges, the PSPC data 
also require a contribution from edges around ~ 0.6 — 0.65 
keV corresponding to edges from O iv-vi. 

Consideration of these maximum allowed ionization 
states for oxygen indicates that the maximum tempera- 
ture of the warm gas is more like T « 5 x 10^ K if the gas 
is isothermal and coUisionally ionized (e.g., Nahar 1999; 
Sutherland & Dopita 1993). The absorption signature of 
this warm gas is not one dominant feature near 0.5 keV 
but is rather a relatively broad trough over energies 0.4 to 
~ 0.8 keV for total optical depths of unity (see Figure 1 
of Krolik & Kallman 1984). 

(We mention that the edge energies we have quoted are 
from Daltabuit & Cox (1972) though Gould & Jung (1991) 
argue that the edge energy for Oi is ~ 10 eV higher. Such 
differences may be relevant for modeling future high reso- 
lution spectra but are unimportant for our present discus- 
sion with the PSPC data.) 

Hence, a proper absorber model needs to consider sev- 
eral edges from different ionization states of oxygen as well 
as edges from ions of C and N. Since the warm gas absorbs 
photons from the hot gas the assumption of coUisional ion- 
ization equilibrium is also not strictly valid nor is it clear 
that the warm gas is fully optica lly thin to its own ra- 
diation as was assumed in §3.6.3 for convenience. Con- 
sequently, the single-edge oxygen absorber that we have 
used throughout this paper (and PAPERl) has been in- 



tended primarily as a phenomenological tool to establish 
the existence and to study the gross properties of the ab- 
sorber which is appropriate for the low spectral resolu- 
tion afforded by the PSPC data. If our basic results are 
confirmed with the substantially higher quality data from 
Chandra and XMM , then it will be appropriate to ex- 
pend the effort to construct rigorous models of the warm 
absorber accounting for many edges and possible radiative 
transfer effects that are not currently available in xspec . 

5.3. Absorber Masses vs Mass Drop-Out 

As suggested in PAPERl this warm ionized absorber 
might be the gas that has dropped out of the putative 
cooling flow during the lifetime of the galaxy or group 
and thus could provide the confirmation of the inhomo- 
geneous cooling flow scenario that has been suggested to 
operate in massive elliptical galaxies, groups, and clusters 
(e.g., Fabian 1994). To make this connection between our 
measurements of oxygen absorption and the cooling flow 
scenario we first estimate the mass of absorbing material 
implied by the measured optical depths. This mass is then 
compared to the cooling flow mass deposition rate inferred 
from ASCA data. 

Before computing the absorber masses some caveats 
must be discussed. First, although the si mple c onstant- 
pressure cooling flow model discussed in §3.6.1 can de- 
scribe the X-ray emission of the hot gas just as well as the 
single-phase model, it does not describe the excess 0.2- 
0.4 keV X-rays presumabl y aris ing from the emission of 
the warm absorbing gas (§ |3.6.3| ). Although at this time 
we cannot exclude the possibility that with a more rigor- 
ous treatment of the absorption and emission properties of 
the warm gas (see end of previous section) that the simple 
cooling flow model could be compatible with the PSPC 
data between 0.2-0.4 keV, it is quite possible that an im- 
portant modifi cati on of the simple cooling flow model is 
required (see §5.5). As a convenient benchmark for com- 
parison to most previous studies we shall consider here the 
mass deposition rates predicted by the constant-pressure 
cooling flow models. 

Second, the oxygen edge optical depths predicted by 
the single-edge models without including the emission 
from the warm gas certainly are overestimates. Recall 
that when including the warm gas emission in NGC 5044 
that T reduces to a value near the 95% lower limit of 
the re sult o btained without including the warm gas emis- 
sion (§ [3. 6. 3 ). Of equal importance, when including addi- 
tional edges the inferred optical depth for each edge de- 
creases, and we expect several edges from different ioniza- 
tion states of oxygen, carbon, and nitrogen to contribute 
(§ ^.2| ). Since the absorption of an edge is not linear in the 
energy of the edge (i.e., Ao{E) = exp[—To{E / Eo)~^] and 
Ai{E)A2iE) ^ Ai+2{E) if the edge energies Ei ^ E2), 
by spreading multiple edges over a large energy range one 
can produce the observed absorption with smaller total 
optical depth than can be achieved with a single edge. 
Con sequently, the single-edge optical depths obtained in 
§3.4 should be considered upper limits. 

Let us now estimate the amount of absorbing mate- 
rial implied by these absorption measurements with these 
caveats in mind, and in particular with the assumption 
that the inferred absorber masses are most likely over- 
estimates. Assuming the optical depths refer to the Oi 
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Table 2 

Absorber Masses and Implied Accumulation Timescales 



Central Bin Total 



Name 


M 

(M^yr-i) 


R 
(kpc) 


(IQiOMo) 


^acc 

(lOiO yr) 


A4ot 
(IQiOMo) 


R 
(kpc) 


Mabs 

(lOiOMf^) 


^acc 

(lOlO yr) 


Mhot 
(10iOM(^) 


N507 

N1399 

N4472 

N4649 

N5044 


18.50+^3:^0 

1 62+"-2° 

2 ■54+0-34 

^•^"-0.60 


20.0 
5.2 
5.2 
5.2 
11.0 


n ai+o-io 
0.81_(, g2 

U.OO_(j 26 

l.69tl;l^ 


QO+O-''^ 
^•■J^_0.23 

12+°-2'' 

16+0-32 
U.lD_o.l3 

04+°-°* 


0-430lj;;g2i 

019+°-°i2 

Q QO-, +0.002 

026+°-°°* 
266+°-°52 


200.0 
90.3 
77.2 
36.7 
143.5 


125 67+599-92 

Q 41+23.11 
^■^^-7.10 
47 26+39-53 

3.70^fgr 

96 c;9+38.17 


g 7Q+40.94 
"•'»_5.29 
c oi +16.78 
&-«i_4.54 
IS fil+21.39 

1.68lfgf 
64+^-1'' 


43.68+12^31 

5 t;ci+0.51 

o fjci+0.36 

n KC+0.19 
0.56_o25 

22.08^4.13 



Note. — Total cooling flow mass deposition rates (Af) and 90% confidence limits inferred from ASCA data (see text). The quoted errors 
on Mabs refer to 95% confidence limits on r ("Fix" in Table |^) and assume solar abundances. The accumulation timescale, tacc = M^i,g/M, 
reflects both the 95% uncertainties on Mabs and the 90% errors on M. The mass of hot gas, A4^t: 'S-'id the 95% errors computed within the 
3D radius are also given. These values of Mabs probably are over-estimates (see beginning of [jp.3[). 



edge, then the measured values of r imply a hydrogen col- 
umn density (assuming a = 5.5 x 10"^^ cm^ at threshold) 
and thus a mass within a projected radius, R, 



A/abs(< R) = 

(7.8 X 10^)(t) 



R 



lOkpc 



0/H 



8.51 X 10-4 



M0, (1) 



where r is the optical depth of the O i edge and O /H is 
the oxygen abundance of the absorber. Although the pro- 
jected mass is larger than the mass within the 3D radius 
r — R, the value of r in equation p) slightly underesti- 
mates the 3D value as discussed in (§p^); i.e. these projec- 
tion effects approximately cancel. The metallicity of the 
hot gas in the central bins for the objects in our sample are 
larger than solar (PAPER2), and thus we expect the same 
for the absorber. Since the oxygen abundance is uncertain 
we shall quote results assuming O is solar and recognize 
that Mabs could be overestimated by a factor of 2-3 in the 
central bin. The expected contribution from carbon and 
nitrogen (§H) to r also reduces Mabs by another ~ 40%. 

In Tableg we give Mabs for NGC 507, 1399, 4472, 4649, 
and 5044 in both the central bin {R = 1') and the total 
mass interior to the largest bin investigated; the edge op- 
tical depths used refer to the single-phase models sii ice th e 
cooling flow models give entirely consistent values (§ 3.6.1 ). 
The mass deposition rates, M, listed in the second column 
are determined from the accumulated ASCA data within 
radii of r ~ 3' — 5'. (The ASCA spectra place much tighter 
constraints on the total M than do the ROSAT spectral 
data. ) The results for NGC 1399, 4472, and 5044 are taken 
from |Buote (19991 ). For NG C 507 and 4649 we re-a nalyzed 
the data sets as prepared in Buote fc Fabian (1998 ) andfrt- 
ted cooling flow models analogously to that done in Buott 



(1999D. That is, the spectra were fitted with (1) a constant 
pressure cooling flow component, (2) an isothermal com- 
ponent representing the ambient gas, and (3) for NGC 
4649 an extra high-temperature bremsstrahlung compo- 
nent. Since the cooling flow model assumes constant pres- 
sure and neglects the gravitational work done on the cool- 
ing gas, the value of M is an upper limit. This over- 
estimate is typically ^ 30% (e.g., the agreement of differ- 
ent cooling flow models in Allen et al 2000b). 



Let us focus on Afabs and the accumulation time, iacc = 
Afabs/M, within the central bin (i? = 1') where the mea- 
sured optical depths are most significant and the fits are 
most clearly improved when the edge is added; i.e. we 
consider the results for the central bin to be most reliable. 
Assuming an age of the universe, tage = 1.3 x 10^°yr, ex- 
amination of Table ^ reveals that within the central bin 
^acc ^ (0.1 — 0.5)iage usiug the best-fitting values or a 95% 
upper hmit of iacc ~ (0.3 - 0.6)iago; NGC 5044 actually 
has smaller values, but if the second bin is included (which 
has intrinsic absorption just as significant as the inner bin) 
then iacc is consistent with the values quoted above. These 
accumulation timescales are a sizeable fraction of iagc, and 
thus A/abs within the central bin(s) can account for most, 
if not all, of the mass deposited by the cooling flow over 
the lifetime of the flow; the exact value depends on pre- 
cisely when the cooling flow begins and whether M varies 
with time. 

The total absorbing masses have large errors within the 
95% confidence limits. Although the best-fitting values for 
iacc are typically larger than tago, the 95% lower limits are 
^ tago for all but NGC 4472. We reiterate that we consider 
these values at the largest radius to be less secure than the 
central bin(s) because the fits do not clearly require the ad- 
dition of the oxygen edge outside the inner 1-2 bins. Nev- 
ertheless, the result for NGC 4472 is striking and deserves 
comment. Clearly the approximation of a spherically sym- 
metric, relaxed cooling flow is invalid for i? ^ 3' because 
the isophotal distortions sug gest a strong interaction with 
the surrounding Virgo gas (Irwin & Sarazin, 1996) and 
thus the estimate of M unlikely applies at larger radii. If 
the large value of Mabs at large radius is confirmed then 
another mechanism must have produced the warm gas in 
NGC 4472. 

Hence, within the central 1-2 bins (i? ^ 10 — 20 kpc) 
where the model constraints are most secure we conclude 
that the absorbing mass inferred from the oxygen edges 
can explain most (and perhaps all) of the mass deposited 
by a cooling flow over the age of the system. If the edges 
also apply at larger radii (as we have assumed), then all of 
the deposited mass (except for NGC 4472) can be easily 
explained by the inferred absorbing mass. These quali- 
tative conclusions still apply if the oxygen edge optical 
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depths are really closer to their 95% lower limits as dis- 
cussed near the top of this section. 

We mention that systems without strong cooling flows 
will not have had sufficient time to accumulate the ~ 
IO^Mq within i? ~ 1' to produce detectable soft X-ray 
emission. Thus, the "Very Soft Components" found in 
galaxies with low ratios of X-ray to optical luminosity un- 
likely arise from warm gas deposited in a cooling flow and 
instead probably reflect the collective emission from X-ray 
binaries (e.g., Irwin & Bregman 1999). 



5.4. Constraints from the Optical and FUV 

Since coUisionally ionized gas at temperatures of 10^^^ 
K emits many strong lines at optical and ultraviolet wave- 
lengths (e.g., Pistinner & Sarazin 1994; Voit & Donahue 
1995), we consider whether the large amounts of absorbing 
material implied by the intrinsic X-ray absorption (Table 
violate published constraints on line emission in the 
optical and UV spectral regions. The best published con- 
straints available in the optical are for Ha from studies 
of extended ionized gas in the centers of elliptical galaxies 
(e.g., Trinchieri & di Alighieri 1991; Goudfrooij et al 1994; 
Macchetto et al 1996). In most cases the emission line 
gas is only detected within r ^ 20" which is significantly 
smaller than the central 1' used in our analysis. 

The object in our sample where Ha has been det ected 
out to the larges t angular radius is NGC 5044 
chetto et al (1996) measure F(Ha) = 1.4 x 10"-^^ erg cm 



i 



Mac 



^ within R = 0.5'. We can estimate the temperature at 
which the Ha emission implied by Mabs (Table |^) within 
R = 1.0' equals the observed fiux. We take the predicted 
Ha line intensity at peak temperature from Pistinner & 
Sarazin and the temperature dependence of ~ 7^-2-5 fj-Qj^j 
inspection of Figure 7 of Voit & Donahue. After account- 
ing for the different region sizes we find that the required 
temperature is w 1.5 x 10^ K using the best-fitting Afabs, 
although when using the 95% lower limit on Afabs we find 
that r « 0.5 X 10^ K. 

Similar results hold for NGC 1399, 4472, and 4649 
although the comparison is less certain because of the 
larger aperture corrections. If no aperture correction is 
made for NGC 5044 then the implied temperatures rise 
to T w 2.5 X 10^ K at best fit and T « 1.7 x 10^ K at 
the 95% lower limit. If we consider also that A/abs in Ta- 
ble 1^ is over-estimated because the oxygen abundance is 
larger than solar 1.5Zq - see PAPER2) and carbon 
and nitrogen contribute ^ 40% to the measured optical 
depths (|5.3|), we obtain T « 1.8 x 10^ K at best fit and 
T « 1.2 X 10^ K at the 95% lower fimit (again without 
aperture correction). Therefore, the published constraints 
on Ha are satisfied if T ^ 2.0 x 10^ K. 

Stronger lines from warm gas are expected to appear 
in the UV, but Hopkins Ultraviolet Telescope (HUT) ob- 
servations detected no significant emission lines in NGC 
1399, 4472, and 4649 (Ferguson et al 1991; Brown, Fergu- 
son, & Davidsen 1995). It is unfortunate that the strongest 
emission lines for temperatures T ^ (2 — 3) x 10^ K are 
Ov(1218A), Ovi(1034A), and Nv(1240A) which appear to 
be lost in the background geocoronal emission (e.g.. Figure 

1 of Ferguson et al 1991, though see below for O VI). How- 
ever, the lines Civ(1549A), Oiv(140lA), and Neiv(1602A) 
are also strong and uncontaminated by geocoronal emis- 



To determine the gas temperature at which, e.g., the 
Oiv(140lA) flux would not violate the published UV con- 
straints we estimate that the O iv flux would have to be less 
than 10% of the continuum considering the error bars 
on the spectrum of NGC 1399 (Figure 2 of Ferguson et al 
1991). This fimit corresponds to a flux of ^ 1.4 x 10^^^ 
erg cm~^ . We take the predicted Oiv line intensity 
at peak temperature from Pistinner & Sarazin and assume 
the temperature gradient above the peak falls similarly to 
that displayed for C iv in Figure 7 of Voit & Donahue. Us- 
ing the best-fitting Afabs (and accounting for the smaller 
HUT aperture) we find that a temperature of at least 
3 X 10^ K is required, though the 95% l ower limits on 
Afabs (which are probably more reafistic - § 3.6.3 ) coupled 
with the oxygen abundance and C /N issues as above indi- 
cate the limit is more conservatively ~ 2 x 10^ K. Similar 
limits are obtained for the other lines and for the HUT 



spectra of NGC 4472 and 4649 ( [Brown et al, 1995| ). 

Our procedure of requiring the line fluxes to be less than 
10% of the continuum may re sult in limits that are too 
restrictive. Dixon et al (1996) have estimated 2a upper 
limits on the Ovi(1034A) intensity from M87 which has a 
HUT spectrum very similar to NGC 1399. Their 2a upper 
limit on the Ovi fiux within a 1' circle is ^ 1 x 10~^° erg 
cm~^ s~^ . If the warm gas has a temperature of 3.2 x 10^ 
K corresponding to the peak temperature for O vi, then 
the observed limit is comparable to the Ovi emission ex- 
pected from the warm gas of NGC 1399 when using the 
95% lower limit on Afabs- Hence, if the HUT results for 
M87 apply to NGC 1399 (as they appear to), then the 
predicted Ovi emission agrees with the limits, especially 
if the temperature is not precisely at the peak temperature 
for Ovi. 

5.5. Theoretical Issues 

Although the hypothesis of warm, mostly coUisionally 
ionized, gas apparently can explain the X-ray observations 
and the matter deposited by the cooling flows, this model 
has serious theoretical difficulties which must be overcome 
before it can be considered a viable model: 

1- Twarm < ^viriai- The temperature of the warm gas 
(^ 0.1 keV) is less than the virial temperatures of 
the halos (~ 1 keV) implying the gas is not ther- 
mally supported. How does the warm gas support 
itself in the gravitational fields of these systems? 

2. t'^^^^^'^ <^ ^coof''''- The cooling time of the warm 
gas is very short, even more so than the hot gas in 
the central regions. How are large quantities of gas 
maintained at these temperatures? 

Most likely these problems can only be solved with a 
substantial modification of the standard cooling flow sce- 
nario. Clearly an additional energy source and very pos- 
sibly the important role of magnetic fields will have to be 
considered. A small number of models have been proposed 
throughout the years which consider these issues, though 
they are not well developed and do not at present make 
detailed predictions regarding the issues (1) and (2) above. 
These models have generally been designed to completely 
suppress cooling flows, and thus will have to be modified to 
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allow some cooling of the hot gas. We now briefly review 
some of the candidate models. 



Bin hey (1996) and ICiotti fc Ostriker (20001) have pro- 



there exist estir nates using minimu m energy argu ments 
for NGC 1399 (|Killecn et al, 1988D, NGC 4472 (|Ekers| 



fc Kotanyi, 1977D , and NGC 4649 ( [Stanger fc WarwickJ 
1986) which give consistent results: S J> 50 — 100/iG at 

10)nG at r 



posed feedback from the central black hole as a promising 
means of inhibiting the cooling of hot gas in the central re- 
gions of galactic cooling flows. In this model whenever the 
black hole accretes a sufficient amount of gas to stimulate 
nuclear activity, the accompanying radiation stimulated 
by the accretion heats up the hot gas and prevents further 
cooling. This is supposed to be a cyclical process such 
that the AGN phase is sufficiently rare to be consistent 
with the lack of nuclear activity in most cooling flows. If 
the AGN feedback energy does not completely suppress 



coolinig of the hot gay but iiiHtead mei'cly pi'evenLH moyt 
the gas fiuni cuuling down Lu LcnipeiaLuies beluw 10'''"'^ 
K, this might be able to explain issues (1) and (2). We 
also note that most of the cooling gas is not ex pected to 
reach the central black hole in standard models ( [Brighcntj 
fc Mathews, 19991 ). 

Another energy source proposed to exist in the centers 
of cooling flows is that from the reconnection of tangled 
magnetic field lines. It has been known for some time 
that small seed magnetic fields can be amplified within 
a cooling fiow to produce sizeable fields at the centers 
(e.g., Soker & Sarazin 1990; Lesch fc Bender 19 90; Moss 
fc Shuku rov 1996; Mathews fc Brighenti 1997). ^oabi et 
al (1998) show that reconnection energy can significantly 
reduce the cooling rates and may account for si gnificant 
amounts of warm (T = 10^-s K) gas. Similarly, f?ormarj 
fc Meiksen (1996| ) propose a two-phase model to recycle 
warm and hot gas along magnetic flux loops also resulting 
in much lower rates of mass deposition. 

Both the AGN and magnetic field reconnection models 
would explain the warm gas phase as a non-equilibrium 
configuration where the warm gas is continuously diluted 
by heating processes and replenished by cooling from the 
hot phase. If we consider the (unlikely) possibility that 
the warm gas is a long-lived equilibrium phase (see, e.g., 
Fabian 1996) then magnetic pressure is probably the most 
viable non-thermal process which can support the gas. 
Mathews fc Brighenti (1999| ) describe a possible equilib- 
rium model of the warm gas as the the outer envelopes of 
low mass stars forming in a cooling flow. 

The details of h ow the magnetic fie ld would support the 
gas are uncertain. Daines et al (1994) suggest that the cool 
gas blobs would be anchored to the hot gas by the mag- 
netic fields, and thus the pressure support would actually 
come from the hot gas. However, inspection of Table 
reveals that Mabs ~ (1 — 10)Mhot within the central bins 
indicating that the hot gas could not support the cool gas. 
(The 95% lower limits on Mabsare considered here.) At 
larger radii it may be possible for the hot gas to support 
the cool gas. 

Whatever the details of the magnetic support the con- 
dition of hydrostatic equilibrium requires that ^ 
eMabsGMgrav?""^- Usiug the values for Mabs within the 
central 1' bin quoted in Table ^ and the gravitating masses 
obtained from previous ROSAT studies (David et al 1995; 
Kim fc Fabbiano 1995; Rangarajan et al 1995; Irwin fc 
Sarazin 1996) we find that B ^ 100/iG is required within 
a 5 kpc radius. 

Estimates of the magnetic field strengths from radio po- 
larization analyses of these galaxies are lacking, although 



the centers and B ^ (5 
B r~.r 



0.5'. Assuming 



(Mathews fc Brighenti , 1997) then these obser 



vations imply {B} ~ 5/iG when averaged over a 1' circle. 
Since the observations only set lower limits the expected 
(B) ~ 100/iG fields are consistent with the observations. 

Interestingly, the n eed for B ^ 100/iG in coolin g flows 
has been suggested by Brighenti fc Mathews (1997) on en- 
tirely different grounds. In their analysis of the g ravitating 



mass distributions of NGC 4472, 4636, and 4649 Brighent 



fc Mathews (1997) find in every case that the gravitat 
ing mass determined from the X-ray analysis falls below 
that estimated from stellar dynamics for r <^ (few)kpc. 
If i? 100/iG within the centers of these systems then 
the X-ray and stellar dynamical masses agree. (A similar 
result holds for NGC 1399 as well - W. Mathews 2000, 
private communication.) 

Finally, we mention it may be useful to consider a model 
for the warm gas that does not originate from cooling out 
of the hot phase. One such possibility is offered by the ma- 
terial continuously ejected by the stars. This material is 
injected into the ISM with low energy and is shock-heated 
up to the virial temperature of the halo. Although the- 
oretical argumen ts suggest that the heating is very rapid 
(Mathews, 1990), more detailed calculations are required 
to rule out the possibility that in fact the transition is grad- 
ual and could give rise to an observable phase of warm gas 
consistent with the X-ray observations. 

6. CONCLUSIONS 

From deprojection analysis of the ROSAT PSPC data 
of 10 cooling flow galaxies and groups with low Galac- 
tic columns we have detected oxygen absorption at the 
2a /3a level intrinsic to the central ^ 1' in half of the sam- 
ple: NGC 507, 1399, 4472, 4649, and 5044. The data for 
the other systems are insufficient to place interesting con- 
straints on the absorption profile but are consistent with 
substantial absorption. We modeled the oxygen absorp- 
tion as a single edge (rest frame E — 0.532 keV) which 
produces the necessary absorption in both the PSPC and 
ASCA data for E ^ 0.5 keV without violating the PSPC 
constraints over 0.2 ~ 0.4 keV for which no significant 
excess absorption is indicated. Assuming the absorber is 
collisionally ionized gas we infer a temperature of 10^~^ K 
from consideration of the possible edge energies consistent 
with the PSPC data. 

The intrinsic oxygen absorption reconciles the long- 
standing problem of why negligible column densities for a 
foreground absorber with solar abundances were inferred 
from ROSAT data whereas large columns were obtained 
from ASCA and other instruments with bandpasses above 
^ 0.5 keV. Moreover, since the absorption is confined to 
energies above ~ 0.5 keV there is no need for large columns 
of cold H which are known to be very inconsistent with the 
negligible atomic and molecular H measured in galactic 
and cluster cooling flows (e.g., Bregman et al 1992; O'Dea 
et al 1994). 

In most of the galaxies and groups we have found that 
single-phase and cooling flow models cannot explain the 
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X-ray emission in the soft ( 0.2-0. 4 keV) energy channels 
of the ROSAT PSPC data (§ ^.6.3| ). That is, when TVh of 
the standard absorber model is freely fitted it is found that 
Nil < N^^^ in the central bins of most systems, with NGC 
4472 and NGC 5044 showing the most significant soft ex- 
cesses. If we model this soft emission as coronal gas we 
obtai n temperatures 10^^^ K in excellent agreement with 



those inferred from the energy ranges of the absorption 



edges. 

Hence, the sub-Galactic column densities are consistent 
with a direct detection of the emission from the intrinsic 
absorbing gas. The agreement between the temperatures 
inferred from the emission and absorbing properties of the 
warm gas lends strong support to the ionized gas model. 
In contrast, dust can not explain the excess soft X-ray 
emissi on. (Other problems exist with the dust hypothesis 
- see jj5j|.) 

Our simple estimates of the amount of absorbing matter 
implied by our single-edge absorption measurements are 
consistent with the total amount of matter expected to 
have been deposited by a cooling flow ( |5.3| ). With the ar- 
rival of higher quality data from Chandra and XMM more 
accurate estimates should be made which account for a 
range of absorbing edges and possible radiative transfer 
effects in the warm gas. 

We have examined the theoretical difficulties associated 
with attributing the absorption to warm ionized gas and 
have discussed some candidate models that may be able to 
account for these problems. Future detailed calculations 



are required to assess the viability of these models (§5.5). 

Fortunately on the observational front it will be very 
easy to verify the intrinsic oxygen absorption with new 
Chandra and XMM data. The XMM (EPIC) and Chandra 
(ACIS-S) CCDs both extend down to 0.1-0.2 keV and have 
substantially better energy resolution than the PSPC. Ob- 
servations with these instruments can easily test our pre- 
diction for warm gas in both absorption and emission. The 
grating spectrometers of XMM and Chandra have even 
better energy resolution (but smaller effective area) and, 
in principle, might detect individual edges. 

We emphasize that the absorption signature of the warm 
gas is expected to be a relatively broad trough over ener- 



gies 0.4 ~ 0.8 keV, and thus future Chandra and XMM 
observations will not see a single sharp feature. The most 
straightforward means to confirm our results will be to 
reproduce the sen sitiv ity of A^h to i^min for a standard 
absorber model (§3.3). To obtain the properties of the 
absorber (e.g., temperature and abundances) a model for 
the soft X-ray o pacity such as that described by Krolik & 
Kallman (1984| ) must be compared to the new data. In 



so doing the emissi on from the warm gas must also be ac- 
counted for (§ |3.6.3 ), and thus it is very important that the 
detector bandpass extend down to ~ 0.1 keV which it does 
for the Chandra and XMM CCDs. 

As discussed in § ^.4| optical and FUV constraints imply 
a lower limit of T ~ 2 x 10^ K. It is possible that precise 
measurements of Oiii(5007A) could refine this limit, but 
since this line peaks at T ~ 0.8 x 10^ K its emissivity is 
already falling rapidly at T ^ 2 x 10^ K. Future high res- 
olution FUV spectroscopy of the O VI (1034A) line (peak 
temperature 3.2 x 10^ K) with, e.g., FUSE may be able to 
place additional interesting constraints on the warm gas if 
its temperature does not exceed T « 5 x 10^ K. 

It should be remembered that to infer the properties of 
the warm gas from X-ray observations the absorption and 
emission spectrum arising from warm gas must be disen- 
tangled from Galactic absorption and the emission from 
hot plasma. Since (if confirmed) the warm gas almost cer- 
tainly represents the mass deposited by an inhomogeneous 
cooling flow (§^), the hot gas at each radius should also 
emit over a range of temperatures. Hence, the thermo- 
dynamic state of the X-ray emitting plasma appears to be 
very complex in the central regions of the (X-ray) brightest 
galaxies and groups, and the analogous results for A1795 
presented in PAPERl suggest the same applies for galaxy 
clusters. 
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